Continuous quantum measurement with particular output: pure wavefunction 

evolution instead of decoherence 



Alexander N. Korotkov 

GPEC, Departement de Physique, Faculté des Sciences de Luminy, Université de la Méditerranée, 13288 Marseille, France 
and Nuclear Physics Institute, Moscow State University, Moscow 119899, Rússia 

(February 9, 2008) 



00 ■ 

■ 

o\ : 

3 ■ 

<: 

r> : 



> : 
\o . 

(N ■ 

o ■ 
oo : 
o . 

00 ■ 

-i— > ■ 
G ■ 
cd 

3 : 



We consider a continuous measurement of a two-level sys- 
tem (double-dot) by weakly coupled detector (tunnel point 
contact nearby). While usual treatment leads to the grad- 
ual system decoherence due to the measurement, we show 
that the knowledge of the measurement result can restore the 
pure wavefunction at any time (this can be experimentally 
verified). The formalism allows to write a simple Langevin 
equation for the random evolution of the system density ma- 
trix which is reflected and caused by the stochastic detector 
output. Gradual wavefunction "collapse" and quantum Zeno 
effect are naturally described by the equation. 

PACS numbers: 73.23.-b, 03.65.Bz 



The problem of quantum measurements has a long his- 
tory, however, it still attracts considerable attention and 
causes discussions and even some controversy, mainly 
about the wavefunction "collapse" postulate (see, e.g., 
Ref. |ï|). Among different modern approaches to this 
problem let us mention the idea to replace the collapse 
postulate by the gradual decoherence of the density ma- 
trix due to the interaction with the detector j|] and the 
approach of a stochastic evolution of the wavefunction 
(see, e.g., ||-||])- The renewed interest to the measure- 
ment problem is justified by the development of the ex- 
perimental technique, which allows more and more ex- 
perimental studies of the quantum measurement in òp- 
tics and mesoscopic structures. The problem has also the 
close connection to the rapidly growing fields of quantum 
cryptography and quantum computing (see, e.g., 0). 

In the recent experiment |?J with "which-path" in- 
terferometer the suppression of Aharonov-Bohm inter- 
ference due to the detection of which path an electron 
chooses, was observed. The weakly coupled quantum 
point contact was used as a detector. The interference 
suppression in this experiment can be quantitatively de- 
scribed by the decoherence due to the measurement pro- 
cess §p. 

We will consider somewhat different setup: two quan- 
tum dots occupied by one electron and a weakly coupled 
detector (point contact nearby) measuring the position 
of the electron. The decoherence of the double-dot den- 
sity matrix due to the measurement has been analyzed 
for this setup in Refs. [UQ. In the present letter we 
answer the following qüestions: how the detector current 
looks like (as a function of time) and what is the proper 



double-dot density matrix for particular detector output. 
We show that the models of point contact considered in 
Refs. PHlCfl describe an ideal detector. In this case the 
density matrix decoherence is just a consequence of ig- 
noring the measurement result. The observer who follows 
the detector output knows the pure wavefunction at each 
moment of time. Moreover, a "mixed" density matrix can 
be gradually purified during the measurement. 

Similar to Ref. let us describe the double-dot sys- 
tem and the measuring point contact by the Hamiltonian 
H = Hdd + T-Ípc + Umt , where Udd = {e/2)(c\ci - 
c 2 c z) + H(c\c2 + c\c\) is the Hamiltonian of the double- 
dot, Hpc = J2l E l a Í a l+12r E raÍa r +J2l, r T ( a i a l+ a l a r) 

describes the tunneling through the point contact {H and 

T are real), and TLi n t = J2i r c 2 c ^( a i a i + a i a r)i i-e. 
the tunneling matrix element for the point contact is T 
or T + AT depending on which dot is occupied. So, 
the average current 1\ = 2irT 2 pip r e 2 V/h flows through 
the detector when the electron is in the first dot (V is 
voltage), while the current is I2 = 1 1 + AI = 2tt(T + 
AT) 2 pip r e 2 V /h when the second dot is occupied. 

We make an important assumption of weak coupling 
between the double-dot and the detector (the better term 
would be "weakly responding" detector), 



|A/| < J = (h + J 2 )/2, 



(1) 



so that many electrons, N > (I /AI) 2 , should pass 
through the point contact before the observer is able 
to distinguish which dot is occupied. This assumption 
allows the classical description of the detector, namely 
the coherence between the quantum states with different 
number of electrons passed through the detector can be 



neglected 12 



The decoherence rate r<j = (y/ï/e — ^fh/e) 2 /2 of 
the double-dot density matrix a(t) due to the measure- 
ment by point contact has been obtained in Ref. ||. 
In the weakly-coupled limit (0) it can be replaced by 
Td = (A/) 2 /8eJo or by expression 



r d = (AI) 2 /4S Iy 



(2) 



where Si = 2eio is the usual Schottky formula for the 
shot noise spectral density Si. Equation (Q) has been also 
obtained in Refs. ]9|-[ïï|] for the quantum point contact as 
a detector, the difference in that case is <Sj = 2e/o(l — T) 
where T is the transpareney of the channel flÏ3| (while 
above we implicitly assumed T< 1 [H)- Notice that 
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the decoherence rate (g) was derived in Refs. |8|-|ll|] with- 
out any account of the information provided by the de- 
tector, implicitly assuming that the measurement result 
is just ignored. Now let us study how this additional 
information affects the double-dot density matrix. 

We start with the completely classical case when there 
is no tunneling between dots (H — 0) and the initial 
density matrix of the system does not have nondiagonal 
elements, 012(0) = 012 (i) = 0. We can assume that the 
electron is actually located in one of the dots, but we just 
do not know exactly in which one, and that is why we use 
probabilities crii(O) and 022(0) = 1—0-11(0). The detector 
output is the fluctuating current I(t). The fluctuations 
grow when we examine I(t) at smaller time scales, so 
some averaging in time ( "low-pass filtering" ) is necessary. 
Let us always work at sufficiently low freqüències, / ~ 
r- 1 « Si/e 2 . 

The probability P to have a particular value for the 
current averaged over time r, (7) = L I(t)dt, is given by 
the distribution 

P((I),r) =an(0)Pi((7),r) +o 22 (0) P 2 «7),r), (3) 
Pi((I),r) = {2^DY 1 ' 2 exp [-((/) - hf /2D] , (4) 

where D = Sj/2r. After the measurement during time 
t the observer acquires additional knowledge about the 
system and should change the probabilities au according 
to the Standard Bayes formula for a posteriori probability. 
Hence, 

ou (r)= ou (0)exp[-((7) -h) 2 /2D] 
x {<7n(0)exp[-((7) -h) 2 /2D] 

+ o 22 (0) exp[-((7) -I 2 ) 2 /2D]y\ 
o 2 2(t)= 1 - on(r). (5) 

Now let us assume that the initial state is fully co- 
herent, 012(0) = a/oh (0)022(0) (while still 77 = e = 0). 
Since the detector is sensitive only to the position of elec- 
tron, the detector current will behave exactly the same 
way. So, after the measurement during time r we should 
assign the same vàlues for au (r) as in Eq. (g). Then the 
upper bound for the nondiagonal element is 



Re 012 (r) < |oi 2 (t)| < v /oii(t)o 2 2(t). 



(6) 



If the actual measurement result is disregarded, then 
the upper bound for 012 can be calculated using the prob- 
ability distribution of different outcomes given by Eq. (||) 
and the upper bound (||) for each realization: 



(Reoi 2 (r)) < / ^a n (r)a 2 2(T)P({I),T) d(I) 

= ^11(0)022(0) cxp[-(A/) 2 r/4ft]. (7) 

This upper bound exactly coincides with the result given 
by decoherence approach (|^). This fact forces us to ac- 
cept somewhat surprising statement that Eq. (pt) gives 



not only the upper bound, but the true value of the 
nondiagonal matrix element, i.e. the pure state remains 
pure after the measurement (no decoherence oceurs) if 
we know the measurement result. 

Simultancously, we proved that the point contact de- 
tector considered in Refs. [js|-[Ï0t causes the slowest possi- 
ble decoherence of the measured system, and hence rep- 
resents an ideal detector in this sense. In contrast, the 
result of Ref. [|Ï5) shows that a single-electron transis- 
tor with large tunncl rcsistances and biased by relatively 
large voltage is not an ideal detector (the non-ideal de- 
tector has been also considered in Ref. pH). 

If the initial state of the double-dot is not purely co- 
herent, |o"i2 (0) < ^/ ÍJ ii(0)o22(0), it can be treated as 
the statistical combination of purely coherent and purely 
incoherent states with the same on(0) and 022(0), then 

Oi 2 (r) = 012(0) [a 11 (T)a 22 (T)/a 11 (0)a 22 (0)} 1 / 2 . (8) 

Eq. (||) together with Eq. (gj) is the central result of the 
present letter. 

Equations (|)-(g) and (§) 

can be used to simulate the 
detector output I(t) and the corresponding evolution of 
the density matrix. For example, in the Monte-Carlo 
method we should first choose the timestep r satisfying 
inequalities e 2 /Si <C r <C Si /(AI) 2 and draw a random 
number for (7) according to the distribution (||). Then 
we update on(í) and a 22 (t) using this value of (7) and 
repeat the procedure many times (the distribution for the 
current averaged over the interval Aí = r is new every 
timestep because of changing au (t) which are used in Eq. 
(^)). The nondiagonal matrix element can be calculated 
at any time using Eq. (||). 

This Monte-Carlo procedure is equivalent to the fol- 
lowing nonlinear Langevin-type equation for the density 
matrix evolution (equation for on is sufHcient): 



Ò11 =1Z, 1Z 



2AI ,r, \ 
-O11O22 —pr~ L·'W 



= -O11O22 



2A7 



Si 

022 - on 



A7 + Ç(t) 



(9) 
(10) 



where the random process £(£) has zero average and the 
low frequeney spectral density Sç = Sj. The second ex- 
pression for 1Z allows to simulate the measurement while 
the first one can be used to calculate the density matrix 
for given 7(t) (that can be done easier using Eq. (||)). No- 
tice that Eq. (M) is closely connected with the Quantum 
State Diffusion approach of Refs. ^g]. 

Figure [ï] shows a particular result of the Monte-Carlo 
simulation for the symmetric initial state, On(0) = 
022(0) = 1/2. Thick line shows the random evolution of 
on(í). Equation ( |Ï0| ) describes the gradual localization 
in one of the dots (first dot in case of Fig. |ï|). Let us de- 
fine the typical localization time as t; oc = 2Si/(AI) 2 (we 
choose the exponential factor at on = 022 = 1/2). Then 
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it is exactly equal to the time Tdis = 2S//(A/) 2 neces- 
sary for the observer to distinguish between two states 
(defined as the relative shift of two Gaussians (0) by two 
Standard deviations), and t; oc = where T d = T^ 1 . It 
is easy to prové that the probability of final localization 
in the first dot is equal to cxii(O), because au(r) averaged 
over realizations is conserved. 

The detector current I(t) basically follows the evolu- 
tion of <7u{t) but also contains the noise which depends 
on the bandwidth. The dashed line in Fig. shows the 
current averaged over the "running window" with the du- 
ration Ai = Si/(AI) 2 while the thin sòlid line is current 
(/} averaged starting from t = 0. Notice that our result 
for I(t) directly contradicts the point of view presented 
in Ref. Q. 

Now let us consider the general case of the double- 
dot with non-zero tunneling H. If the frequency Í2 of 
"infernal" oscillations is sufficiently low, Í7 = (AH 2 + 
e 2 ) 1 / 2 /h <C Si/e 2 , we can use the same formalism just 
adding the slow evolution due to finite H (the product 
Çlrioc is arbitrary). The particular realization can be ei- 
ther simulated by Monte-Carlo procedure similar to that 
outlined above [now update of a\2 (t) using Eq. (||) should 
be necessarily done at each timestep] or equivalently de- 
scribed by the coupled Langevin equations 



àn = -à 22 = {-2H/h)lm(a 12 )+TZ, (11) 
H 

fi " 2<7n<722 



it ÍH (T22 — fil 
àl2 = T"0"12 + -r-(o"n — CT22) + — 7i(Tl2 



-7rf (J 12, 



(12) 



where 7<j = for an ideal detector (see below). 

Figure |^ shows particular results of the Monte-Carlo 
simulations for the double-dot with e = H and dif- 
ferent strength of the interaction with an ideal detec- 
tor. The electron is initially located in the first dot, 
cii(O) = 1. The dashed line shows the evolution of <7n 
with no detector. Notice that because of e ^ 0, the ini- 
tial asymmetry of the electron location remains in this 
case for infinite time. When the interaction with detec- 
tor, C — h(AI) 2 1 SjH, is relatively small (top sòlid line), 
the evolution of en is close to that without the detec- 
tor. However, the electron gradually "forgets" the initial 
asymmetry and the evolution can be described as the 
slow variation of the phase and amplitude of oscillations 
(recali that the wavefunction remains pure). 

When the coupling with the detector increases, the 
evolution significantly changes (middle and bottom 
curves in Fig. ||). First, the transition between dots 
slows down (Quantum Zeno effect). Second, while the 
frequency of transitions decreases with increasing inter- 
action with detector (at sufficiently strong coupling), the 
time of a transition decreases, so eventually we can say 
about uncorrelated "quantum jumps" between states. 

In a regime of small coupling with detector, C « 1, 
the detector output is too noisy to follow the evolution 



of au and, correspondingly, only slightly affects the oscil- 
lations. On contrary, when C > 1 the detector accurately 
informs about the position of electron and simultaneously 
destroys the oscillations. 

Equations can be generalized for a 11011- 

ideal detector, T d > (A/) 2 /4S*/ (as in Refs. 0|]§), 
which gives an observer less information than possible 
in principle. Let us model it as two ideal detectors "in 
parallel" with unaccessible output of the second detec- 
tor. Thcn the information loss can be represented by 
the extra decoherence term — 7d<7i2 in Eq. ( |Ï2| ) where 
ld = — (AI) 2 /4Si. The limiting case of a nonideal 
detector is the detector with no output (just an environ- 
ment, AI = 0) or with disregarded output. Then the 
evolution equations (|ïl·|)-(|Ï2]) reduce to the Standard de- 
coherence approach. 

For nonideal detector it is meaningful to keep our 
old definition of the localization time, r; oc = túís — 
2S7/(AI) 2 while r d < 2n oc . So, we consider localization 
time not as a real physical quantity (that is meaningless 
because observer cannot check it) but as a quantity re- 
lated to observer's information. Similarly, the effective 
decoherence time is defined as r' d — 7^ 1 . 

Equations ©-© with the term TZ given by Eq. @ 
can be used to obtain the evolution of the density ma- 
trix in an experiment provided the known detector out- 
put I(t) and initial condition (7jj(0). Notice that even if 
the initial state is completely random, o\\ = 021 = 1/2, 
(T12 = 0, the nondiagonal matrix element appears during 
the measurement due to aequired information, so that 
sufficiently long observation with an ideal detector leads 
to almost pure wavefunction. Such a purification of the 
density matrix described by Eqs. (|Tl·|)— (|Ï2[) is analogous 
to the localization at H = 0. 

Let us briefly discuss the philosophical aspect of the de- 
veloped formalism. The statement that the pure wave- 
function remains pure during the continuous measure- 
ment by an ideal detector (with known result) may seem 
surprising at first, however, we easily recognize the direct 
analogy to the "orthodox" situation of a "sharp" mea- 
surement (the wavefunction is pure after the "collapse"). 
Anothcr important point is that the density matrix is 
in some sense observer-dependent. For example, if two 
observers have different level of access to the detector 
information (as, for example, in the model of nonideal 
detector considered above), then the density matrix for 
them will be different. Nevertheless, the observer with 
less information still can use his density matrix for all 
purposes, just his predictions will be less aceurate. So, 
instead of the "actual" density matrix, it is better to dis- 
cuss only the "accessible" density matrix, that is fully 
consistent with the "orthodox" (Copenhagen) point of 
view. 

If the knowledge of the detector output is not used 
in the experiment, the decoherence approach is suitable. 
On contrary, one can devise an experiment in which the 
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subsequent system evolution depends on the preceding 
measurement result; then the proper description is given 
by Eqs. 

For example, let us consider the double-dot with H = 
and fully coherent symmetric initial state. According 
to our formalism, after the measurement by an ideal 
detector during time r (most interesting case is t < 
t; oc ) the wavefunction remains pure but becomes asym- 
metric (Eqs. (||) and (||)). To prové this, for exam- 
ple, an experimentalist should switch off the detector 
at t = t, reduce the barrier between the dots (create 
finite H), and create the energy difference e = [(1 — 
4|cri2 1 2 ) 1/2 - l]HReai 2 /\cri2\ 2 ; then after the time pe- 
riod Ai = [tt — arcsin(Im<Ti2 híl/ H)]/fl the electron will 
be moved to the first dot with the probability equal to 
unity, that can be checked by the detector switched on 
again. Alternatively, using the knowledge of 0"íj(t) an 
experimentalist can exactly prepare the ground state of 
the coupled double-dot system and check it, for example, 
by the photon absorption. Another experimental idea is 



FIG. 1. Thick line: particular Monte-Carlo realiza- 
tion of un evolution in time during the measurement of 
uncoupled dots, H = 0. Initial state is symmetric, 
fií (0) = 022(0) = 1/2, while the measurement leads to grad- 
ual localization. Initially pure wavefunction remains pure at 
any time t. Thin line shows the corresponding detector cur- 
rent (7) averaged over the whole time interval starting from 
t = while the dashed line is the current averaged over the 
running window with duration Si/(AI) 2 . 

FIG. 2. Random evolution of au (particular Monte-Carlo 
realizations) for asymmetric double-dot, e = H, with the elec- 
tron initially in the first dot, <7n(0) = 1, for different strength 
of coupling with detector: C = h(AI) 2 /SiH =0.3, 3, and 30 
from top to bottom. Dashed line represents C — (unmea- 
sured double-dot) . Increasing coupling with detector destroys 
the quantum oscillations (while wavefunction remains pure at 
any i), slows down the transitions between states (Quantum 
Zeno effect), and for C ^> 1 leads to uncorrelated jumps be- 
tween well localized states. 



to start whh cumplelely laiidum sLalc uf llic duublc-dut 
witll íiniLe II and Lhcii giadually (inusL iiiLeiesLing case 

is ÍÍT/ oc < 1) obtain almost pure wavefunction using the 
detector output I(t) and Eqs. ®-((ï|). The hnal test 
of the wavefunction is similar to that considered above. 

An experiment of this kind can verify the formalism 
developed in the present letter. While such an experi- 
ment is still a challenge for the present-day technology, 
we hope that it can be realized in the near future. 

In conclusion, we developed a simple formalism for the 
evolution of double-dot density matrix with account of 
the result of the continuous measurement by weakly cou- 
pled (weakly responding) point contact. The formalism 
is suitable for any two-level system measured by weakly 
coupled detector. 
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